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ABSTRACT 



We have constructed color-magnitude diagrams (CMDs) for ten globular clusters in 
the halo of the nearby spiral galaxy M33 based on Hubble Space Telescope Wide-Field 
Planetary Camera 2 observations in the F555W (~V) and F814W (~I) filters. These 
data reveal the morphology of the HB and allow us to estimate the cluster metallicity 
using the shape and color of the red giant branch. The principal result we report 
herein is that 8 of the 10 clusters possess exclusively red HB morphologies yet their 
metallicities are as metal-poor as [Fe/H] = —1.6. Indeed, these 8 clusters present 
basically only giant branch clumps reminiscent of intermediate age star clusters in 
the Magellanic Clouds. In addition, two of the clusters form a second parameter pair 
which have similar metallicities but very disparate HB types. Under the assumption 
that cluster age is the global second parameter, the average age of halo globular 
clusters in M33 appears to be a few Gyr younger than halo clusters in the Milky Way. 
Using the observed properties of HB stars in M31 and M33 along with published 
main-sequence turnoff ages for the globular clusters in the Milky Way, LMC, SMC, 
and the Sagittarius dwarf spheroidal (Sgr), we attempt to sketch the early formation 
history of these galaxies. This indicates that the Milky Way, M31, M33, the LMC, 
and Sgr all experienced their first epoch of cluster formation soon after the Big Bang. 
Three to four Gyr later, the SMC began to form its first generation of clusters; the bulk 
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of the M33 clusters formed later still. We note that the halo clusters in M33 formed 
over a much larger time period than those in the Milky Way and M31. 



Subject headings: galaxies: star clusters, halos, formation 



1. Introduction 



The globular clusters of the Milky Way (MW) probe the earliest formation epochs of our 
Galaxy. Identifying similar clusters in other galaxies allows us to study the properties of these 
galaxies soon after their formation. In addition to the MW, a number of other Local Group 
galaxies contain globular clusters; among these are the Sagittarius dwarf spheroidal galaxy (Sgr), 
the Large and Small Magellanic Clouds (LMC and SMC), the Fornax dwarf spheroidal galaxy, 
M31, and M33, which is the topic of the present work. We have constructed CMDs for a sample 
of halo globular clusters in M33 in order to estimate their metallicities and HB morphologies. By 
comparing these properties with those of globular clusters in other Local Group galaxies, we hope 
to construct a formation chronology for the Local Group. 



2. Sample Selection, Observations and Reduction 

We are interested in the oldest star clusters in the halo of M33. As a result, our target 
selection was primarily based on the integrated color and kinematics of the clusters. In particular, 
we inspected Fig. 2 of Schommer et al. (|1991| ), which shows the difference between the cluster 
velocity and the disk velocity as a function of the integrated cluster color. This difference is 
essentially zero for the blue clusters but begins to show a significant dispersion redward of 
[B — V) ~ 0.6, indicating that clusters redder than this value are likely to be located in the halo 
of M33. We were also interested in sampling as much of the spatial and radial extent of M33 as 
possible. All of these points lead us to our target list of 10 M33 clusters whose properties are 
tabulated in Table 1. Assuming (m — M)q = 24.5 for M33 (van den Bergh 1991), these clusters 



have -9<-My<— 7, placing them squarely among the Galactic halo globular clusters in luminosity 



(Harris 1997). In addition, our sample includes all clusters deemed by Christian & Schommer 



11983|) to have ages > 4 Gyr. Previous studies, including Christian & Schommer (|1983[) , Cohen et 



al. (|1984 ), Christian & Schommer flL98fl ), Schommer et al. (|1991|) , and Brodie & Huchra (|1991|) , 
indicate that, while M33 possesses clusters with a wide age range, the clusters we selected are 
among the oldest clusters in this galaxy. To sum up, our target list of 10 M33 clusters are similar 
to Galactic halo globular clusters in terms of their luminosities, colors, kinematics, measured 
metallicities, and implied ages as they were known at the time. 



The observations used in the present work were obtained with the Wide-Field Planetary 
Camera 2 on HST as part of our cycle 5 program (GO-5914). Each cluster was centered in the 
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PCI chip. We observed each target for 4 orbits, yielding a total exposure time of 4800 s in F555W 
and 5200 s in F814W. 

After reprocessing using the latest 'recommended' calibration exposures, the program PCI 
frames were reduced using the procedure outlined by Silbermann et al. ( 1996| ), which describes the 



techniques used by the HST 'Cepheid Distance Scale' key project. Further details are provided in 
Sarajedini et al. ( |1998j ). In summary, empirical point spread functions (PSFs), kindly provided by 
P. Stetson, were fitted to the detected stellar profiles on each CCD frame using the ALLFRAME 
package. The resultant photometry was edited and matched to form instrumental colors. Aperture 
corrections derived from our frames and the photometric transformations of Silbermann et al. 



( 1996 ) were applied to bring the photometry onto the standard system. Along with the typical 
error in the aperture corrections of ~0.08 mag, the estimated uncertainty of ~0.05 mag in the 
transformations leads to a total uncertainty of ~ ±0.09 mag in the photometric zeropoints. 



3. Results 

A preliminary visual inspection of the cluster CMDs revealed a rather surprising result. 
Although we had selected the M33 clusters to be analogous to Galactic globular clusters in every 
respect, all but two of the M33 halo clusters possess extremely red HBs, or red clumps lying 
essentially on top of the cluster RGBs. Color-magnitude diagrams of M33-M9 and M33-C20 
are shown in in Figs, la and lb, respectively. These include stars between 1 and 2.8 arcsec of 
the cluster centers. The solid line is the fiducial sequence of the Galactic globular cluster M5 
(Sandquist et al. 1996] ) which has a metallicity between that of M9 and C20. In addition, the lower 



panels of Fig. 1 show the CMDs of M3 (Buonanno et al. 1994| ), which has a metallicity similar to 



that of M33-M9, and NGC 362 (Harris |1982j ), which is similar to M33-C20 in metallicity. Both of 
these Galactic globulars are classified as 'younger halo' clusters by Zinn (|i"993|) . Fi gure 1 is an 



illustration of the fact that, while two of our M33 clusters (M9 and U77) have HB morphologies 
similar to the 'younger-halo' MW clusters, the majority of the M33 halo clusters (8 of the 
10 in our sample) possess HBs that are significantly redder than any Milky Way 
globular clusters at their metallicity. This is the primary result of the present work. 

In order to quantify the metallicities, we rely on the slope of the RGB. Mighell et al. ( |1998| ) 
have shown that metallicities derived from the RGB slope are minimally affected by the age of 
the cluster. The measurement of the slopes and the transformation to metallicities on the Zinn 



& West (1984) system are both described in Sarajedini et al. (|1998|) . Table 1 lists the measured 
metallicities for our program clusters. Because of the extreme differential reddening in the region 
of M33-R14, its abundance has been visually estimated via a direct comparison with the RGBs of 



Da Costa & Armandroff ( 1990| ), Therefore, the estimated metallicity of R14 is rather uncertain. 



Table 1 also lists the cluster metallicities from Christian & Schommer ( |1988| ; [Fe/H]cs) and 
Brodie & Huchra ( pM| ; [Fe/H] bh), both of which have been derived using integrated spectra. 



Our values are in general intermediate between these other studies, but the differences are not 
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statistically significant. The reddening for each cluster (Table 1) follows from Eq. 1 of Sarajedini 
(1994). The errors in the reddenings do not include the uncertainty in the photometric zeropoint. 



The significance of these values is discussed in Sarajedini et al. (1998). 

To quantify the HB morphologies, we rely upon generalized color histograms of the HB stars. 
For the red HB clusters where all of the HB stars lie in a clump on or near the RGB, we adopt the 
peak of this histogram, i.e. ~(V — I) g , as the mean HB color. For the two clusters with significant 
blue HB populations, we perform an 'intensity-weighted' average of the colors of the red and blue 
peaks. We therefore have quantitative measures of each cluster's metallicity and intrinsic mean 
HB color (i.e. HB morphology). 



4. Discussion 

To facilitate the comparison of these colors with those of Galactic globular clusters and 
keeping in mind that the absolute errors in our photometry are fairly significant, we will examine 
the difference in color between the HB and RGB (ds-v = (B — V) g — (B — V)hb)- Although it 
makes very little practical difference, for the sake of consistency, we have converted the (V — /) 
values of the M33 clusters to {B — V) using a combination of the relation published by Zinn & 
Barnes ( p-996j ) and one described by Sarajedini et al. (1998). The filled symbols in Fig. 2 are the 



observed values of [Fe/H] and ds-v for Galactic globular clusters from Buonanno et al. ( 1997 ) 
and Sarajedini et al. ( |1995| ). In the former case, we have adopted the relation between metallicity 
and dereddened RGB color from Sarajedini & Layden ( [1997 ) to convert {{B — V) Q: hb) to ds-v- 



The plus signs in Fig. 2 are 7 'old' globular clusters in the SMC from the work of Mighell et 
al. (|1998| ). The open circles show the locations of the 9 M33 clusters for which we have reliable 
measurements. The values for M9 and U77 (clusters with a blue component to their HBs) are 
indicated as lower limits because artificial star experiments (described in Sarajedini et al. |1998j ) 
indicate that our observations/reductions are not able to detect an extended blue component to 
the HBs of these clusters, if in reality they possess such a feature. In addition, some of the ds-v 
values for the red HB clusters are slightly negative; this is simply due to measurement errors, 
which could be as much as ~0.05 mag. The dashed lines in Fig. 2 represent the location of 
scaled-solar abundance Zero Age Horizontal Branches for masses of 0.62, 0.66, 0.74, and 0.90 M & 
from the work of Dorman (1992), again converted to ds-v y i a the relation published by Sarajedini 



& Layden ( |1997|) 



Figure 2 shows that the general trend is for more metal-poor clusters to have bluer HB colors, 
and this is corroborated by the theoretical models. However, there is a significant dispersion 
in ds-v at a given metallicity, especially at intermediate abundances. This is the so-called 
"second parameter effect" wherein another parameter, in addition to metallicity, is influencing the 
morphology and thus the color of the HB. Figure 2 suggests that, whatever this second parameter 
actually is, the majority of halo globular clusters in M33 are significantly more affected by it than 
those in the MW or the SMC. In fact, the M33 clusters U49 and M9 are a second-parameter pair 
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(similar to M3 and M13) in which both clusters have similar metallicities, but vastly different 
HB morphologies. Among the many possible candidates for the second parameter, Sarajedini et 
al. ( 1997| ) assert that cluster age is the most likely one (although see Stetson et al. 1996| for an 
alternative viewpoint). If we adopt this assertion, then the M33 globular clusters appear to be 
significantly younger, in the mean, than similar clusters in the MW. 

There is at least one additional piece of evidence that supports the possibility of young age 
for these M33 globular clusters. Several of them appear to have a small number of stars located 
brighter than the first ascent RGB tip. Figure 3 illustrates the CMDs of C20, H38, H10, and R12. 
The photometry has been shifted in color using the reddenings given in Table 1 and in magnitude 
by requiring that (V — I)o,-3 (the dereddened color at Mi = —3) be equal to that predicted 



by Equation 1 of Da Costa & Armandroff (1990) given the metallicity in Table 1. This is done 
instead of utilizing the RR Lyrae magnitude because these clusters have red HB clumps, which 
are generally not at the RR Lyrae level. The solid lines in Fig. 3 are the standard RGB sequences 
of Da Costa & Armandroff ( |199C| ). The presence of stars above the magnitude of the first ascent 
red giant branch tip (M/~-4.0) for clusters with [Fe/H] ^ — 1.0 supports the assertion that some 



of the M33 halo globular clusters are of young to intermediate age (see Mould & Da Costa 1988). 



We note that there were previous indications that perhaps some of these clusters may not 
be as old as globular clusters in our Galaxy (Cohen et al. |1984| ; Searle et al. |198C| ; Christian & 



Schommer |1983|) . To place a lower limit on the ages of these clusters, we have compared the CMD 



of C20 with the theoretical isochrones of Bertelli et al. ( 1994 ). This comparison indicates that the 
ages of the youngest M33 clusters must be older than ~2 Gyr; otherwise the main-sequence turnoffs 
would be bright enough to be detected in our CMDs. Given this fact, we use the nomenclature of 
Hodge ( 1988|) to remind ourselves that "globular" clusters should be characterized by their mass 
(> 10 4 Mq) and not by their age. Estimating an upper limit to the ages of these clusters is less 
straightforward. Certainly, we can say that the two clusters with blue HB components (M9 and 
U77) are likely to be similar in age to the 'younger-halo' MW globulars (such as M3). However, 
an upper limit to the ages of the remaining eight M33 globular clusters, which have red clumps 
located essentially on top of their RGBs, is more difficult to ascertain. From these observations, 



we can conclude that the halo kinematics reported by Schommer et al ( 1991 ) for this sample 
of objects implies that some intermediate- aged M33 clusters are NOT in the main disk of that 
galaxy, which has implications for cluster formation in a wide range of galaxies. We also note in 
passing that, since only 20% of our M33 clusters are bona fide old globular clusters, the specific 
frequency of such clusters in M33, which Schommer et al. (1991) and Bothun ( |1992[ ) have claimed 
is anomalously high (~60-70), may be reduced by a factor of ~5. 



5. Implications for the Formation of the Local Group 



By studying their main-sequence turnoffs, it has become clear that the oldest globular clusters 
in the LMC are similar in age, within ~1 Gyr, to those in the MW (Olsen et al. 1998| ). In 
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addition, Mighell et al. ( [1998] ) have used main-sequence turnoff photometry of the oldest globular 
clusters in the SMC and find that these are 3-4 Gyr younger than the oldest MW clusters. The 
old clusters in the Sagittarius dwarf galaxy (Sgr), M54 and Terzan 8, have been studied by Layden 
& Sarajedini ( |1997| ) and Montegriffo et al. ( |1998] ), respectively. Both of these clusters appear to 
be as old as MW globular clusters at their metallicity. 

We can correlate the findings based on main-sequence turnoff photometry with the observed 
HB morphologies of the oldest clusters in each galaxy discussed above. For example, the oldest 
globular clusters in the MW, LMC, and Sgr all possess HBs that are populated predominantly 
blueward of the RR Lyrae instability strip. In contrast, the oldest clusters in the SMC, which are 
3-4 Gyr younger than those in the MW, LMC, and Sgr, have predominantly red HB morphologies. 
This trend of redder HB morphology with younger age is, in our opinion, simply the second 
parameter effect noted above. 

A number of studies have exploited the superior resolution of HST to construct CMDs for 
M31 globulars (Ajhar et al. [199§ Fusi Pecci et al. [l99|; Holland et al. fll997D ). The collective 



results of these investigations reveal the existence of metal-poor blue HB clusters, which seems to 
indicate that the oldest globular clusters in M31 are probably similar in age to those in the MW 



(see also Djorgovski et al. 1997). Along the same lines, in the halo of M33, we have uncovered two 



intermediate-to-metal-poor clusters with blue+red HBs. Thus, based upon the appearance of Fig. 



2 and the implications of theoretical HB models (e.g. see Lee et al. 1994 ), we can be reasonably 
certain that, even though the MW, M31, M33, LMC, and Sgr are spread across a region of space 
several hundred kiloparsecs in size, they all formed their first generation of clusters at around the 
same time soon after the Big Bang (see also Harris et al. 1997). Ages determined from cluster 
main-sequence turnoffs indicate that 3 to 4 Gyr later, the first globular in the SMC (NGC 121) 
formed. Lastly, the bulk of the M33 clusters (those with purely red HBs) followed several Gyrs 
later still. If we consider only the galaxies with 'halo' populations, we find that the age range 
among the halo clusters in M33 is somewhere between a factor of 2 to 3 larger than those of the 
MW and M31. 
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Fig. 1. — The top two panels show CMDs for M9 and C20, which are two halo globular clusters 
in M33. The solid lines are the fiducial sequence of the Galactic globular cluster M5 (Sandquist et 
al. 1996) offset to match the V(HB) and (V - I) g of M9 and C20. The lower panels show CMDs 
of the Galactic globular clusters M3, which is similar to M33-M9 in metal abundance, and NGC 
362, which is similar to M33-C20 in metal abundance. The M33 CMDs are V vs. V-I while those 
for the Galactic globular clusters are V vs. B-V. 

Fig. 2. — The difference in color between the HB and RGB {ds-v) f° r MW globular clusters (filled 
circles), SMC globular clusters (plus signs), and the M33 halo globular clusters (open circles). The 
dashed lines are the expected locations of zero-age horizontal branch models for 0.62, 0.66, 0.74, 
and 0.90 M Q (from right to left) from the work of Dorman (1992). 

Fig. 3. — Color-magnitude diagrams for four of the M33 clusters in our sample. The solid lines 
are the standard red giant branches of M15 ([Fe/H] = -2.2), NGC 1851 {[Fe/H] = -1.3), and 47 
Tuc {[Fe/H] = -0.7) from Da Costa & Armandroff (1990). 
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Table 1. M33 Globular Cluster Characteristics 



Cluster 


V 


B-V 


i? a 


[Fe/H] cs b 


[Fe/H} BH ^ 


[Fe/H] CM D d 


E(V - I) 


HB Morphology 


U49 


16.25 


0.68 


8.4 


-0.8 ±0.3 


-1.70 ±0.53 


-1.64 ±0.20 


0.07 ±0.02 


Red 


R12 


16.38 


1.03 


4.1 


-1.2 ±0.3 




-1.19 ± 0.24 


0.05 ± 0.03 


Red 


R14 


16.48 


0.98 


2.5 


-1.5 ±0.3 




-1.0 ±0.5 




Red 


M9 


17.12 


0.69 


8.9 


-1.7±0.3 




-1.64 ±0.28 


0.04 ± 0.03 


Red±Blue 


U77 


17.19 


0.67 


5.4 




-1.77 ±0.77 


-1.56 ± 0.30 


0.08 ± 0.03 


Rcd+Bluc 


H38 


17.25 


0.73 


11.1 


-1.5 ±0.3 




-1.10 ± 0.10 


0.04 ± 0.01 


Red 


C20 


17.67 


0.77 


15.3 


-2.2 ±0.3 


-1.25 ±0.79 


-1.25 ± 0.22 


0.03 ± 0.03 


Red 


C38 


18.10 


0.73 


17.5 


-1.2 ±0.3 




-0.65 ± 0.16 


0.04 ± 0.02 


Red 


H10 


18.23 


0.96 


10.5 




-0.91 ±0.90 


-1.44 ±0.26 


0.25 ± 0.03 


Red 


U137 


18.30 


0.83 


12.8 




-0.12 ±0.38 


-0.98 ±0.16 


0.09 ± 0.03 


Red 



a Radial distance in arcminutes from the center of M33. 
b Christian & Schommer (1988) 
c Brodic & Huchra (1991) 

d Metallicity based on the color-magnitude diagram. 
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Figure 1 - Sarajedini et al. 
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Figure 2 - Sarajedini et al. 
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Figure 3 - Sarajedini et al. 



